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ABSTRACT 

We present a detailed calculation of model atmospheres for DA white dwarfs. Our 
atmosphere code solves the atmosphere structure in local thermodynamic equilibrium 
with a standard partial linearization technique, which takes into account the energy 
transfer by radiation and convection. This code incorporates recent improved and 
extended data base of collision induced absorption by molecular hydrogen. We analyse 
the thermodynamic structure and emergent flux of atmospheres in a range 2500 < 
Teff < 60000 K and 6.5 < \ogg < 9.0. Bolometric correction and colour indices are 
provided for a subsample of the model grid. Comparison of the colours is made with 
published observational material and results of other recent model calculations. 

Motivated by the increasing interest on helium core white dwarfs, we analyse 
the photometric characteristics of these stars during their cooling, using evolutionary 
models recently available. Effective temperatures, surface gravities, masses and ages 
have been determined for some helium core white dwarf candidates, and their possible 
binary nature is briefly discussed. 
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1 INTRODUCTION 

White dwarf (WD) stars represent the most common end- 
product of the stellar evolution. On average, WDs have a 
mass of O.6M0, a central density of 10® gr cm~^ and a core 
composed of carbon and oxygen resulting from helium burn- 
ing during the preceding phases of the evolution (D'Antona 
& Mazzitelli 1990). Due to the high degree of electron de- 
generacy in the interior of WDs, the outer nondegenerate 
layers, and specially their atmospheres, control the loss of 
energy into outer space and therefore, the cooling of these 
objects. From a technical view point, detailed atmosphere 
models provide an adequate surface boundary condition for 
the computation of evolutionary stellar models. Besides, the 
atmosphere models are necessary to determine the emergent 
flux from the star, which is required to interpret the lumi- 
nosities and colours of observed WDs. 

Current interest in WDs is motivated by two reasons. 
These stars are candidates for the Galactic dark matter and 
can also be used to estimate the age of the star system they 
belong to. The recent detection of gravitational microlens- 
ing events in the direction of the Magellanic Glouds (Alcock 
et al. 1997, Alcock et al. 1999, Renault et al. 1997) suggests 
that very low mass stars, WDs among others, might consti- 
tute a significant fraction of the mass of the Galatic halo. 
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The luminosity function necessary to study this interpreta- 
tion requires the calculation of cooling tracks using accurate 
atmosphere models. 

The spectral evolution of a WD during its cooling can 
be used to calculate the age of the star. The first detailed 
calculation of WD evolution was performed by Iben & Tu- 
tukov (1984), for the cooling of a carbon-oxygen WD from 
the planetary nucleus stage through the stage of complete in- 
ternal crystallization. Recent works (Hansen 1999; Chabrier 
1999) show that the time scales involved, particularly in an 
advanced stage of the cooling process, depend strongly on 
the atmosphere composition (among other factors such as 
the stellar mass), which through the opacity regulates the 
energy loss from the stellar interior. Thus, for instance, the 
strong opacity from hydrogen molecules causes a slow cool- 
ing of the WDs with hydrogen-rich atmosphere, resulting 
these objects potentially observable for times comparable to 
the Hubble time (Richer et al. 2000). The technique based 
on the observation of a low-luminosity cutoff has been ap- 
plied to estimate the age of the local Galatic disk (Winget 
et al. 1987; Wood 1992) and the age of clusters (von Hippel, 
Gilmore & Jones 1995; von Hippel & Gilmore 2000). 

The majority of WDs have pure hydrogen atmospheres 
due to gravitational settling, which is efficient at removing 
helium and heavier elements from atmosphere toward inner 
layers (Fontaine & Michaud 1979; Muchmore 1984; Iben & 
MacDonald 1985; Althaus & Benvenuto 2000). Convective 
mixing could lead to a helium-rich atmosphere in cool WD 
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when the mass of the superficial hydrogen layer is small 
enough (Koester 1976; Vauclair & Reisse 1977; D'Antona 
& Mazzitelli 1979). In particular, Bergeron, Ruiz & Leggett 
(1997) have reported evidence for convective mixing in cool 
WD atmospheres. However, hydrogen-rich atmospheres are 
observationally most abundant than those of helium or other 
with metals, in part due to the long cooling times. Further- 
more, these atmospheres present a structure which is insen- 
sitive to small admixtures of helium (Hansen 1999) . 

In this context, recent model atmosphere calculations 
(Hansen 1998, 1999; Saumon & Jacobson 1999) suggest that 
WDs of very low luminosity, and therefore very old, become 
bluer and not redder as they age (and their effective temper- 
ature decreases), resulting their detection easier. This being 
the case, an important number of these objects should be 
observed with the new generation of telescopes. 

Another important question of current interest is the 
evolution of WDs with very low mass having helium cores 
(He WD). Such stars, due to the low stellar mass, would 
be formed from some binary systems. Indeed, numerous bi- 
nary configurations of this type have been observed lately 
(Marsh 1995; Moran, Marsh & Bragaglia 1997; Edmonds et 
al. 1999) . In particular, low mass WDs have been detected in 
binary systems composed of a WD star and a pulsar. A par- 
ticularly interesting binary system with these characteristics 
is the millisecond pulsar J1012-I-5307 and its WD compan- 
ion, because it is possible to determine the properties of the 
pulsar from the analysis of the WD (Van Kerwijk, Bergeron 
& Kulkarni 1996). The study of these binaries would allow 
us to determine an upper limit to the neutron-star mass 
(and in this way, to constrain the still uncertain properties 
of the equation of state of the cool nuclear matter) and on 
its age (with implications for the evolution of its magnetic 
field). Obviously, to obtain accurate results from the study 
of these systems it is necessary to have suitable evolution 
models and, in particular, detailed model atmospheres for 
WDs. 

The primary aim of this paper is to present results of 
model calculations from an independent atmosphere code. 
This will be applied to determine accurate cooling tracks 
of WDs, by providing appropriate outer boundary condi- 
tions to evolution models. Our model atmospheres with pure 
hydrogen compositions cover extensive ranges of effective 
temperature (2500 < T^ff < 60000 K) and surface gravity 
(6.5 < log g < 9.0). Although several sets of models available 
in the literature (cf. Bergeron, Ruiz & Leggett 1997; Saumon 
& Jacobson 1999; Hansen 1999; Chabrier 1999) have already 
investigated different parts of this {Teff,g) region, we offer 
here an overall analysis. We have considered it important 
to extend the study to the very low gravities (log g = 6.5) 
which are found in early stages of the WD formation, and 
also in the cool stages of He WDs with very low mass, which 
have received little or no attention so far. In particular, we 
have also computed model atmospheres following the cool- 
ing tracks of He WDs with hydrogen envelopes and several 
masses. 

The results presented here are based on ideal equations 
of state. The incorporation of density effects in our code 
using the occupation probability formalism of Hummer & 
Mihalas (1988) is in progress. However, nonideal effects can 
be neglected for hydrogen WD models with Te// > 4000 K 
(Bergeron, Saumon & Wesemael 1995a). They increase for 



lower Te// and are significant for Te// < 2500 K (Saumon 
& Jacobson 1999). 

In Section ^ we derive explicit analytical expressions for 
the thermodynamic quantities required in the calculations. 
The method used for the model atmosphere computations 
is presented in Section ^, together with the opacities and 
convective theory adopted. In Section ^ we display the at- 
mosphere structure of WDs obtained from our code. Com- 
parisons with observations and other recent model results, 
through the analysis of colour diagrams, are also presented 
there. Our conclusions are given in Section ^ 



2 THERMODYNAMIC QUANTITIES 

The computation of convective flux of energy requires some 
properties of the gas such as the adiabatic gradient and the 
specific heat. To determine them, we consider a mixture of 
H2, H, H"^, e~, and radiation occuping a volume V, in ther- 
modynamic equilibrium at temperature T. In the present 
calculations, species with low abundances can be ignored 
(Saumon & Chabrier 1992). 

Let A*' be the number of hydrogen nuclei or protons. 
Then, the number conservation and charge neutrality con- 
strains demand that the particle density is 

xN 

'V ' 

with 



(1) 



l + Xl 



X2 



(2) 



where xi = nu+ /n and X2 = riH^/n (with n = N/V) are the 
number concentration of H+ and H2. We use the ideal state 
equation to express the gas pressure. The total pressure, 
including the radiation contribution, is then 

„ xNkT 1 4 

P=-^ + -aT^ (3) 

We adopt a simple approach for the total energy in the 
system 



U 



-xNkT + 



kT + kTA- + 



1 



X2 X2 



X2N + arv , 



(4) 



where 6 = Tv/T. The first term on the right-hand side is the 
translational energy, which assumes the particles to behave 
as classical point particles. The second term gives the poten- 
tial energy of the hydrogen nuclei, which increases with the 
dissociation of H2 molecules and the ionization of H atoms 
according to the corresponding binding energies xi and X2. 
For simplicity, we ignore in this term the excited states of 
these species. The zero point of the energy scale is taken to 
be the ground state of H2. The third term represents the 
rotation- vibration energy of H2, where we use a vibration 
temperature r„ = 6100 K. Finally, the last term takes into 
account the energy in the radiation field. 

Expressions for the adiabatic gradient Va = 
d\nT/d\nP\g, and heat capacity at constant pressure 
Cp = dU/dT\p + P dV/dT\j^ca.n be obtain from differen- 
tiation of Equations (^) and We have followed a similar 
procedure to that of Krishna-Swamy (1961; see also Mihalas 
1965) and, therefore, we only show here the results obtained 
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including molecular species in the gas mixture. The expre- 
sion for the adiabatic gradient is 



ai 



with 



ai = 1 + 47 - 61 



and 



a2 = - + I67 + 
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dlnP„ 



+ f)2 



dX2 
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dlnP„ 



X2 , , dxi 



(5) 
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where 



2 kT 'J X 
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For the heat capacity we obtain 
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The partial derivatives of the concentrations xi and X2, have 
been calculated considering that, except in the regime of 
pressure dissociation, H'*' and H2 do not have simultaneously 
high abundances. Then, using Saha equations for the first 
and third reactions in Table 1, result 
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(11) 
(12) 

(13) 
(14) 

(15) 

(16) 



where the dependence of the partition functions with the 
temperature and pressure has been neglected. 

Since the adiabatic temperature gradient forms the ba- 
sis of the Schwarzschild criterion for convective instability, 
this gradient plays a central role in the computation of at- 
mosphere models. Therefore, we check our expression (|^-0) 
by comparing with results from the equation of state (EOS) 



for hydrogen developed by Saumon & Chabrier (summa- 
rized in Saumon, Chabrier & Van Horn 1995, SCVH). This 
EOS, considered the best available at present for hydrogen, 
is based on the free-energy minimization technique (for de- 
tails, see Graboske, Harwood & Rogers 1969; Fontaine, Gra- 
boske & Van Horn 1977; Hummer & Mihalas 1988) and in- 
cludes a sophisticated treatment of nonideal effects. 

Comparisons between values for the adiabatic gradient 
obtained from SCVH and our EOS are shown in Fig. 1. 
The curves correspond to models of WD atmosphere with 
Teff = 2500, 4000, 6000 and 10000 K, and log 5 = 8. Both 
present and SCVH calculations agree quite well except for 
high pressures of the model T^ff — 2500 K, where non- 
ideal effects can be present. For the hottest two models, dis- 
crepancies at low pressure occur beyond the regime where 
SCVH EOS can be applied (4 < logP(dyn cm-^)< 19). In 
particular, the upper layers of the Teff — 6000 K atmo- 
sphere are dominated by atomic hydrogen (see Section |4.l[ ) , 
and a tendency to a value 0.4 for the temperature gradient 
(corresponding to an ideal monoatomic gas) is reproduced 
by our calculations. A partial molecular formation (up to 
nH2/n — 0.12) takes place as one proceeds inward in this 
atmosphere, yielding a dramatic drop in the value of Va- 

Other familiar features can be identified in this figure. 
For instance, Te// — 2500 K model represents an atmo- 
sphere formed mainly by H2. At low temperatures present 
in its superficial layers, the H2 vibrational modes remain 
frozen and then, the adiabatic gradient takes a value close 
to 2/7 = 0.286, i.e. that due to a fiuid of rigid rotating 
molecules with two degrees of rotational freedom and three 
degrees of translational freedom. The adiabatic gradient de- 
creases toward deep layers where the molecular gas is par- 
tially dissociated. 

The low values reached by Va in the Te// ~ 4000 K 
model are due to a marked sensitivity of the molecular abun- 
dance to changes in temperature and pressure. While the 
behaviour of Va in the hottest model {T^ff = 10000 K) is 
produced by different ionization degree at different depth. 

In conclusion, the relative simplicity of our EOS allows 
us to obtain the analytic expressions for the heat capacity 
and the adiabatic gradient necessary to compute atmosphere 
models. This has the advantage of avoiding potential errors 
due to numerical inaccuracies that can be present in a so- 
phisticated EOS, which requires interpolation from tables 
and numerical evaluation of derivatives. 



3 MODEL ATMOSPHERE 

Cooling of WDs is a very slow process so that the atmo- 
sphere structure responds to quasi-static boundary condi- 
tions. Furthermore, the region where the emergent radiation 
is formed has a thickness much less than the stellar radius. 
These properties allow us to choose a static model atmo- 
sphere with the approximation of plane-parallel layers and 
constant surface gravity throughout the atmosphere. 

On the other hand, because the gas densities are very 
high in the layers where the emergent radiation is formed, 
local thermodynamic equilibrium (LTE) is a good approxi- 
mation for the matter state and we adopt that. 

We take the atmophere gas to be composed of H, H2, 
e~, H~, H^, HJ and H^. The main constituents are H2, 
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Table 1. Reactions considered at the chemical equilibrium. 



Reaction ionization/dissociation energy (eV) 



H - 


. H+ 4 


- e 


13.598 


H- 


^ H 4 


- c^ 


0.755 


Ha - 


2H 




4.478 




^ H 4 


- H+ 


2.651 




Ha 


+ H+ 


4.354 



H, H and e~, depending on the physical conditions. The 
other species are usually present as traces, but they can 
play an important role in the absorption and emission of 
radiation. In fact, the H~ ion can be a dominant source 
of opacity in cool atmospheres, whereas Hj" and Hg" are 
important because their presence affects the concentration 
of H~ and e~ . 

To determine the equilibrium concentrations we con- 
sider the reactions given in Table 1. The values of energy 
are those compiled by Lenzuni & Saumon (1992). We use 
the partition functions calculated by Irwin (1981) for H and 
Ha, by Sauval & Tatum (1984) for HJ, and Neale & Ten- 
nyson (1995) for H^. We have checked our chemical equi- 
librium against results from an independent code (Saumon 
1999; based on an equation of ideal state, except in a sharp 
cuttoff in the H partition function depending on the dis- 
tance to the nearest neighbor). Though some differences (at 
most 12%) appear between both calculations (mainly due to 
different ways of computing the Ha partition function), the 
agreement is satisfactory in general. 

3.1 Opacities 

We have consider the most important opacity sources in 
a hydrogen atmosphere. Bound-free and free-free absorp- 
tion coeficients of H have been taken from Mihalas (1978), 
with Gaunt factors from Menzel & Pekeris (1935). Isotropic 
Thomson scattering is assumed for the opacity of free elec- 
trons. We adopted absorption coeficients of the negative hy- 
drogen ion according to John (1988) for bound- free tran- 
sitions, and analytic fits by Gray (1992) based on Bell & 
Berrington (1987) results for free-free transitions. 

Collisions between Ha molecules yield collision-induced 
absorption (CIA), consisting of rototranslational and rotovi- 
brational bands with peaks located in far-red and infrared 
spectral regions. We use Ha-Ha CIA cross sections calculated 
by Borysow and coworkers, including the most recent results 
(Borysow, Jorgensen & Zheng 1997). For increasingly large 
densities, absorption from collisions involving more than 
two Ha molecules becomes important. Three-body collisions 
are no longer negligible for densities p > 0.02 g cm~"^ and 
they dominate the opacity above p > 0.2 g cm~"^ (Lenzuni 
& Saumon 1992). Ignoring terms of higher order, the to- 
tal frequency-integrated CIA coefficient can be expressed as 
(Lenzuni, Chernoff & Salpeter 1991) 

k = fcap' + fcap^ , (17) 

where fca and fcs account for absorption by binary and 
ternary collisions. We assume the ratio fca/fca — 0.05/p me- 
assured experimentally by Hare & Welsh (1958). 

Finally, cross section of Rayleigh scattering from molec- 
ular hydrogen is taken from Dalgarno & Willians (1966). We 



also consider other opacity sources such as Rayleigh scatter- 
ing from H, bound-free and free-free absorption of HJ, and 
free-free opacity of H^. For all of them, we use coefficient 
expressions from Kurucz (1970). 



3.2 Convective Transport 

Energy transport by convection is treated within the mixing- 
length theory (e.g.. Cox & Giuh 1968; Mihalas 1978) in most 
model atmosphere calculations of WDs. Although this phe- 
nomenological formalism involves free parameters, Bergeron, 
Wesemael & Fontaine (1992) have found that model results 
with Tef f below ~ 8000 K are insensitive to the assumed pa- 
rameterization. In addition, they found that the predicted 
emergent fluxes show dependence on the efficiency of con- 
vection in the range Teff ~ 8000-15000 K. However, it is 
difficult to decide which is the most appropriate parame- 
terization for WDs. Preliminary analyses of the effects of 
convective efficiency on atmospheric parameters of ZZ Ceti 
[^and other lukewarm DA stars, suggest that MLl and MLS 
parameterizations (following the nomenclature of Fontaine, 
Villeneuve & Wilson 1981) could be inadequate (Wesemael 
et al. 1991, Bergeron et al. 1992). In this context, the ML2 
version of the mixing-length theory has been considered a vi- 
able alternative. The results presented here use this version 
of convection. 



3.3 Numerical Method 

The model atmosphere may be calculated using an iteration 
procedure over linearized structure equations in terms of 
pertubations of the temperature and density (or pressure). 
However, we adopt a linearization method for the tempera- 
ture alone (Gustafsson 1971), the density and pressure struc- 
ture being obtained from a given temperature distribution 
by integration of the hydrostatic equilibrium equation. This 
procedure yieds similar results to a complete linearization 
method, but has a much shorter computing time (Saumon 
et al. 1994; SBLHB). Additional economy of execution time 
is obtained from formulating the radiative transfer in terms 
of moment equations, with closure using variable Eddington 
factors (Auer & Mihalas 1970). 

To calculate a model atmosphere, the equations of ra- 
diative transfer and constant flux condition are linearized as 
functions of the temperature. The equations are discretized 
and organized following the Rybicki scheme (for details see 
Gustafsson 1971 and Gustafsson & Nissen 1972). The cal- 
culation of a model starts with a given temperature distri- 
bution from which the hydrostatic equilibrium equation is 
integrated, and the opacities and thermodynamic quantities 
derived. Then, transfer and energy constrain equations are 
solved and we obtain a temperature correction. The whole 
process is iterated to convergence. The Eddington factors 
are evaluated in each step from a Feautrier solution of the 
radiative transfer, using known values of the source function. 



t WDs of DA class (i.e. with Balmer lines and no He I or metals 
present, and therefore, having H-rich atmospheres) showing non- 
radial pulsations. 
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3.4 Model convergence properties 

We remark on some points about the performance of the cal- 
culations. Although the correction procedure adopted here 
has proved effective and stable, certainly it is not as stable as 
when the energy transfer is purely radiative. Numerical diffi- 
culties can appear to solve the coupled equations of transfer 
and conservation energy when both convective and radiative 
modes are present (Gustafsson 1971; Grenfell 1974; Koester 
1980; Bergeron et al. 1991; SBLHB). These difficulties are 
not exclusive to partial or complete linearization methods, 
since they have also been reported by Mihalas (1965) using 
a modification of the temperature-correction method pro- 
posed by Avrett and Krook (1963). 

We found that numerical instability can occur mainly 
in atmospheres of WDs with low Te// and high gravity. 
In these models, the convection becomes extremely efficient 
and, therefore, the temperature gradient is nearly adiabatic 
(V — Va < 10~^ in convective deep layers). Then, the tem- 
perature correction oscillates yielding a large convective flux 
in one iteration, and zero convection in that following for the 
same layer, according to the Schwarzschild instability crite- 
rion. Among the different ways that have been proposed to 
stabilize the iterations (see above references), we found con- 
vinient a similar method to the one suggested by Koester 
(1980). This consists of setting Fc = F — Frad instead of 
Fc = G when V < Va for a layer which would be convective 
in the converged model. 

On the other hand, a good initial temperature distri- 
bution decreases the number of iterations necessary and in- 
creases the range of convergence of the model in the (7^//, 
log g) plane. For purely radiative models, one can start with 
an approximate solution of the gray atmosphere problem, 
e.g. the Eddington distribution (Mihalas 1978). However, 
when the convection is present, this class of distributions 
usually leads to an enormous convective flux which may 
have destabilizing consequences on the temperature correc- 
tion method. A more suitable starting model is obtained 
when in the convective zone where Fc > F, the temperature 
gradient is changed and forced to yield Fc — F (Nordlund 
1974). This procedure yields a suitable initial temperature 
stratiflcation in convective models, from which to apply the 
rigorous temperature correction scheme. 

Finally, as in the report of SBLHB, we found that a dis- 
continuity develops in the temperature distribution of con- 
verged models with Te// « 4400-5500 K at logg ^ 7.5, 
extending to lower effective temperatures when the grav- 
ity decrease. Such as discussed further by Bergeron et al. 
(1995a), this effect is due to competition between opacities 
of H2 molecule and hydrogen negative ion. This convergence 
problem can be solved using an appropriate form of the ra- 
diative equilibrium equation (SBLHB), but unfortunately 
needs a large number of iterations. 



4 RESULTS AND DISCUSSION 

We have calculated a grid of models covering an extensive 
range of effective temperatures 2500 < Tef / < 60000 K and 
surface gravities 6 < logg < 9. In these computations we 
used 76 depth points to solve the hydrostatic equilibrium 
and energy tranfer equations in the range — 6 < logr < 1.5. 



The continuum radiation field was computed at 87 to 107 
frequency points depending on the effective temperature of 
the model. The points were also chosen taking into account 
the location of photometr ic b andpasses used in the colour 
computations (see Section 4.3). 

The temperature-correction procedure was repeated un- 
til the temperature changes were everywhere AT/T < 0.1%, 
for most models. In hot models, T^// ^ 14000 K, which 
are radiative or weakly convective, less than 10 iterations 
were sufficient to reach the convergence, yielding deviations 
from flux constancy less than 0.1 per cent everywhere. For 
lukewarm models, 9000 Si Te// ^ 14000 K, convergence of 
AF/F < 1% is obtained with about 10 iterations. As the 
convection becomes more efficient, Te// ^ 9000 K, the num- 
ber of iterations increases to about 40. In the range of tem- 
peratures where the onset of the H2 formation takes place, 
we need to arise the iteration number to approximately 150, 
due to the change of the boundary conditions in energy bal- 
ance equations (see SBLHB). 

4.1 Atmospheric structure 

Strongly convective flux transport develops in the inner lay- 
ers of cool and lukewarm WD atmospheres. This occurs 
preferentially in zones where hydrogen is partially ionized 
or dissociated, which yields a signiflcant drop in the value 
of the adiabatic temperature gradient favouring the convec- 
tive instability. Fig. 2 shows the fraction of flux carried by 
convection in atmospheres with surface gravities log g — 6.5 
and 8.0. It is also shown for log g — 6.5 models the regions 
where H2, H and H^ in turn dominate the gas. Convection 
is absent from the atmosphere in the hotter early stages of 
a WD. According to our calculations, a superflcial convec- 
tion zone appears as the star cools below Te// « 12000 K 
for WDs with low surface gravities (log g = 6.5) and below 
Te// ^ 16000 K for more massive stars (log g = 8). The 
convective zone extends over the region where part of the 
emergent continuum radiation is formed (0.1 ^ t ^ 1). 

For cooler models, approximately Te// < 10000 K for 
log g = 6.5 and Te// < 11000 K for log g = 8.0, the ef- 
ficiency of convection increases and a considerable fraction 
of the fiux (more than 90 per cent) is carried by convection 
in the innermost layers (r ^ 10). Radiation always domi- 
nates in the optically thin layers t ^ 1. However, thanks 
to high densities and molecule formation, a little convective 
flux persists up to very transparent layers, r ~ 10"'^, in at- 
mospheres around Te// ~ 4000 K. The computed models 
predict a reduction of convection for very cool atmospheres 
Te// ~ 2500 K, but this result may be affected by non-ideal 
contributions to the EOS. Slight pressure dissociation of H2 
can be expected at large depths in these models, affecting 
the profile of the convective fiux as consequence of changes 
in the adiabatic gradient (Saumon & Jacobson 1999). 

For log g = 8.0 models, the convection structure dis- 
played in Fig. 2 is very similar to that showed by Bergeron 
et al. (1991, Fig. 2) in the common region (Te// > 5000 K). 
As a difference, our results show as the molecular formation 
(not included in that work) strongly affects this structure 
below Teff = 6000 K, increasing the convection toward su- 
perflcial layers. 

Results about the run of several physical variables are 
shown in Figs 3-5 for selected models with different effec- 
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tive temperatures at log g = 6.5 and 8.0. For the common 
models at log <? = 8, a good agreement is obtained from our 
T{P) profiles compared to those by Bergeron, Wesemael & 
Fontaine (f 995a, Fig. 2, T^ff = 4000 - f 0000 K) and Saumon 
& Jacobson (f999. Fig. f, Teff = 2500 and 3000 K). The 
only differences are found in the transparent layers of hot 
models (Te// = 8000 - 10000 K), where the line blanketing 
(considered by Bergeron et al.) yields there slight changes 
in the temperature distribution. 

Our sets of models with log g = 6.5 and 8.0 present sim- 
ilar aspects except that higher pressures and densities are 
found in each model atmosphere with increasing g (Figs 3- 
5) . In terms of T^f / , it is possible to distinguish roughly three 
classes of atmospheric structures. One corresponds to hotter 
models, T^// ^ 12000 K at log g = 6.5 and T^// Z 14000 
K at log g = 8.0, which are purely radiative or with a su- 
perficial zone weakly convective. They show rather smooth 
temperature and density gradients without additional signif- 
icant details. Throughout these atmospheres, most material 
is ionized. 

As Te// decreases, a new atmospheric structure takes 
place. The models with T^ff « 4500-12000 K at log g = 6.5 
and Teff « 4500-14000 K at log g = 8.0, cover the range 
from the onset of convection with little flux transported by 
it to the models where the molecular formation becomes 
important. In the present group of models, the relation- 
ship between temperature and pressure exhibits (Fig. 3) 
a clear transition from a very small gradient in radiative 
outer layers to a steep gradient in the top of the convective 
zone. In contrast, the run of density with pressure becomes 
relatively shallow in the transition region between radia- 
tive and convective zones (Fig. 5). In this sense, the model 
{Tef f , log g) — (12000, 6.5) is interesting, because it shows a 
slight inversion of the density in approximately log P — 4.7. 
It can be interpreted as result of a too fast temperature rise 
with increasing P (see Fig. 3), which is compensated by a 
decrease of the density (P oc pT) . 

Finally, the coolest models studied, 2500 < T^ff 4500 
K, show some proper features. They are characterized by 
molecular formation which, as mentioned earlier, increases 
the convective efflciency and extends the convective zone 
toward very transparent regions (Figs 2 and 4), reaching a 
maximum at Teff = 4000 K. This effect has also been found 
in very low mass stars with log g — 5 (SBLHB). At such low 
Teff, our models show a uniform temperature in the upper 
radiative layers and a rather sudden rise with depth in still 
transparent regions (Fig 4). Deep inside the convective zone 
results a low temperature gradient, which is very close to 
the adiabatic value. 

Fig. 5 shows that very cool atmospheres extend to ex- 
tremely high density and pressure regimes (note that the top 
and bottom layers of all models are located at log r = — 6 
and 1.5, respectively). Non- ideal effects would actually be 
very small in models with relatively low gravity log g — 6.5, 
but it can become important in the deep layers (below 
p ~ 10~^) of very cool models with log g — 8.0. Never- 
theless, the T{P) structures of the atmosphere models are 
away from the region where Saumon & Chabrier (1992) pre- 
dict the plasma phase transition of hydrogen (Fig. 3) . 

Compared to hot and lukewarm models analysed above 
(except at T^ff around of the onset of convection), the tem- 
perature distributions of very cool atmospheres {Teff < 



4500 K) are particularly sensitive to the value of the gravity. 
The surface temperature decreases and the gradient in the 
(r, T) plane increases as the gravity is raised (Figs 3 and 



4). This is a consequence of the CIA opacity (Section 3.1) 
which, being a coUisional process, becomes stronger in atmo- 
spheres with high gravity (therefore more dense) increasing 
the cooling of upper layers. 



4.2 Emergent flux distribution 

To illustrate the behaviour of the WD spectra in the 
{Teff,g) plane, we display in Figs 6a-6f the emergent flux 
corresponding to some selected Teff values for log g = 6.5, 
8.0 and 9.0. The importance of different opacity sources 
can be noted in Fig. 7, where individual contributions to 
the total monochromatic opacity are shown for some rele- 
vant physical conditions. These correspond to a Rosseland 
optical depth r = 1 in atmospheres with log g = % and 
Teff = 60000, 10000, 6000 and 2500 K (from 7a to 7d, re- 
spectively). Fig. 7d shows, with heavy lines, the opacity con- 
tributions at r = 0.1 for the 2500 K model. 

Emergent fluxes for the three highest values of Teff 
clearly show the effect of the atomic hydrogen absorption, 
which is a dominant source at such high temperatures (Figs 
6a-c and 7a-b). Among the spectra present in Fig. 6, only 
the Tef f — 60000 K models show an important radiation 
flux emerging from the Lyman continuum. 

At Teff = 15000 K and log <; = 9, an important convec- 
tive zone (with a maximum of 40% of the flux transported 
by convection) develops between log r ~ —1.5 and 1.3. This 
strongly affects the Balmer continuum decreasing Lyman 
and Balmer jumps, with respect to lower gravity models 
(Fig. 6b). This sensitivity offers the possibility of an obser- 
vational test of the convection theory even with relatively 
low-quality data. 

Bound-free and free-free process of H~ and atomic 
Rayleigh scattering start to take place as the effective tem- 
perature diminishes to 10000 K. The most prominent feature 
in the spectra of these atmospheres is the Balmer jump, 
which markedly weakens for large gravities. We note that 
these changes in the Balmer jump can also be reproduced 
by the presence of helium in the atmosphere and an increas- 
ing of its abundance (Wegner & Schulz 1981; Bergeron et al. 
1991). 

A large variety of opacity sources are found in the atmo- 
spheres of lukewarm WDs (Fig. 7c), the H~ opacity being 
the most important. Since H~ bound-free opacity is domi- 
nant and weakly frequency dependent, the emergent flux of 
Tef f — 6000 K model is relatively featureless and resembles 
more a blackbody distribution (Fig. 6d). Furthermore, it is 
relatively insensitive to g. 

Fig. 7d emphasizes the importance of the molecular 
opacity at r = 1 in a (Tef f , log g) — (2500,8.0) model. As 
one proceeds outward in the atmosphere, the opacity be- 
comes completely dominated by CIA process and Rayleigh 
scattering (heavy lines in Fig. 7d, corresponding to r = 0.1), 
due to a marked decrease of the abundance of all species but 
the H2 molecule. Besides, as the temperature decreases, sep- 
arate vibrational bands are clearly identifiable in the CIA 
opacity. Note that the decrease of the density (from r — 1 
to r = 0.1) does not affect the molecular Rayleigh scattering 
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contribution to the total opacity per gram, because the H2 
remains as the most abundant component. 

Synthetic spectra obtained for Te// ~ 2500 and 4000 
K (Figs 6e-f) show a well-known result: the effect of CIA 
opacity is largest for WDs of low T^ff and high gravity, 
corresponding to a large abundance of H2 and high density, 
and, therefore, a large number density of H2-H2 pairs. The 
effect of the four vibrational bands of H2 (Fig. 7d) is clearly 
visible in the spectrum of models with Te// = 2500 K. This 
strong absorption over infrared spectral regions coincides 
with the peak emission of a blackbody at the same Te//- As 
a consequence, the stellar flux is forced to emerge at shorter 
wavelengths, around the minimum in the opacity caused by 
the H2 CIA and the Rayleigh scattering. 



4.3 Colour indices 

From the emergent flux provided by the atmospheric model, 
we can derive magnitudes and colours. We have calculated 
broadband colour indices using the Johnson-Cousins BVRI 
and Johnson-Glass JHK response functions of Bessell (1990) 
and Bessell & Brett (1988), respectively. The magnitude m 
corresponding to a bandpass with tranmission function 
was computed from 



-2.5 log 



47r 



(f) 



dX 







S^dX 



+ c„ 



(18) 



where Hx is the monochromatic Eddington flux from the 
model atmosphere for a WD of radius R and distance D from 
Earth. We adopted the calibration constants Cm given in 
Bergeron, Ruiz & Leggett (1997, BRL), yielding magnitudes 
on the Carnegie image tube system. A set of colour indices is 
obtained from the difference between magnitudes, resulting 
independent of the radius and distance of the star. Fig. 8 
shows the location of the fllter transmitions over two illus- 
trative emergent fluxes. We have also calculated the bolo- 
metric correction at the visual using the expression derived 
by Bergeron, Wesemael & Beauchamp (1995b, Eq.[3]) 

/>oo 

BC = 2.5 log / HxSxdX - 10 log T^ff + 15.6165 . (19) 
Jo 

In view of increasing attention on very low mass WDs 
in the last few years, we report the colours and bolometric 
corrections for a subsample of our grid, corresponding to low 
gravity (log g = 6.5) and cool (Te// = 4000-8000 K) models 
(Table 2). We also show the results for a canonical value 
logs = 8.0. 

Fig. 9 shows the {V — R, R ~ I) two-colour diagram 
for all models calculated. The solid lines connect models 
with same gravity and log g — 6.5 (0.5) 9.0 (from top to 
botton). Also displayed is a subset of 44 DA stars (filled 
circles) and 45 non-DA stars (open circles) observed by BRL 

and having complete BVRI- JHK photometric data ^ By 
comparison, we include the hydrogen atmosphere sequence 
at log 3 = 8 of Bergeron (2000; diamonds), and that of 
Saumon & Jacobson (1999; triangles) covering Te// < 4000 
K. The cooling sequence of Richer et al. (2000; dashed line) 



t The sample contains WDs from the catalog of McCook & Sion 
(1987) reclassified DA by BRL. 



for a WD with carbon-oxygen core, hydrogen atmosphere 
and mass 0.6 Mq is also shown. 

Our models predict that all stars with hydrogen atmo- 
spheres above T^// ~ 4500 K, fall on a line of almost con- 
stant surface gravity with little dispersion. The DA stars 
plotted in the figure remain very close to the synthetic colour 
sequences over the range 3500 ^ T^ff ^ 8500 K. In contrast, 
the non-DA stars reveal more dispersion, although most of 
them seem to follow the hydrogen sequences. Only a few non- 
DA stars (marked on the figure) show colour indices strongly 
affected by other non-hydrogen sources of opacity. Among 
them, the very low luminosity WD LHS 3250 has a DC clas- 
sification (continuous spectrum without lines deeper than 
5%) and its atmospheric composition is not known (Harris 
et al. 1999). BPM 27606 is a DQ star (ie. with carbon fea- 
tures) whose spectrum presents very strong C2 Swan bands 
(BRL). LP 790 - 29 is also a DQ star. 

Fig. 9 shows that the theoretical sequence at log g = 
8.0 of Bergeron (2000) is well reproduced by our models 
(this occurs also for other log g values, not shown here). For 
the coolest atmospheres {T^ff < 3500 K), our sequence at 
log (? = 8 becomes bluer than that of Saumon & Jacobson 
(1999). Nevertheless, the three independent models found an 
equivalent location for the turnoff in R — I. This turnoff is a 
consequence of the onset of the H2 formation (which starts 
the CIA and molecular Rayleigh opacities), and therefore it 
varies with the surface gravity. The largest values oi R — I 
correspond, at instance, to Te// = 3500 K at log g = 6.5 and 
Te// = 4000 K at log g = 9.0. In agreement with Saumon 
& Jacobson, we found that several colour indices show the 
presence of a turnover at low effective temperatures, with 
B — V being an exception (see below). 

In addition, the cooling curve of Richer et al. (2000) 
follows the other sequences of models and the sample ob- 
served for Te/ / > 4000 K, but crosses several sequences with 
different g for lower Te// where a constant or smoothly in- 
creasing gravity should be expected as the star cools. We 
also note a certain irregular behaviour where our grid pre- 
dicts that colours are insensitive to changes in the gravity 
(Te// > 4500 K). 

The [B ~V,V ~ K) two-colour diagram is displayed in 
Fig. 10. There is a close agreement between Bergeron's and 
our models at log g — 8.0, except for the hot models. A 
similar turnover inV — K colour (reported earlier by Mould 
& Liebert 1978) is predicted by both sets of models. Never- 
theless, significant departures of our sequence from that of 
Bergeron can be observed below B — V ~ 0.3 (Te// > 7500 
K) . These discrepancies presumably have their origin in the 
hydrogen line opacity and the occupation probability for- 
malism (Hummer & Mihalas 1988) used by his code and not 
included in our calculations. 

The diagram also reveals a certain difficulty of the the- 
ory in reproducing the observed values in the B — V colour 
at low effective temperatures. In fact, there is an increasing 
difference between observed and theoretical colour indices 
above B — V ~ 0.7. BRL have interpreted this in terms of 
a pseudo continuum opacity at the Lyman edge (not taken 
into account in the atmosphere codes), which affects the flux 
in the B magnitude (for details see their paper). 

Fig. 11 illustrates that the sample of observed stars in 
the {V ~ I ,V ~ K) diagram is well reproduced by our model 
grid, which is also in good agreement with the theoretical 
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Table 2. Broadband colour indices of the log g = 6.5 and 8.0 models 
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+0.648 


-0.029 


-0.076 


-0.398 


4500 


8.0 


+0.898 


+0.581 


+2.160 


+0.583 


+0.271 


+0.044 


-0.403 


5000 


8.0 


+0.767 


+0.499 


+2.027 


+0.500 


+0.341 


+0.105 


-0.298 


5500 


8.0 


+0.644 


+0.423 


+1.709 


+0.423 


+0.300 


+0.085 


-0.198 


6000 


8.0 


+0.538 


+0.360 
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+0.056 


-0.154 


6500 


8.0 


+0.440 
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+1.201 
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+0.028 


-0.152 


7000 


8.0 


+0.362 


+0.260 


+1.000 
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+0.007 
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7500 


8.0 


+0.298 


+0.223 


+0.822 


+0.220 
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8.0 


+0.243 


+0.192 


+0.660 


+0.186 


+0.136 


-0.024 


-0.246 



sequence of Bergeron. It is also observed that a marked sen- 
sitivity of y — ii" and V ~ I colours to the effective tempera- 
ture and surface gravity in models beyond the turnover. The 
reason of this is that / and K bandpasses lie just where the 
H2 CIA has strong contributions from two major vibrational 
transitions, while the V bandpass receives a flux excess that 
is forced to emerge in the optical spectral range. One would 
expect that this sensitivity in the WDs with hydrogen at- 
mospheres could be very useful for cosmochronological pur- 
poses. 



4.4 White dwarfs with low mass 

Helium core white dwarf (He WD) stars are considered the 
final product of the evolution of some close binary systems. 
The mass of these objects should be smaller than that re- 
quired for degenerate helium ignition, ~ 0.5 AIq (Mazzitelli 
1989), thus, resulting very low mass WDs. They have been 
detected in large surveys (Bragaglia et al. 1990; Bergeron, 
Saffer & Liebert 1992; Bragaglia, Renzini & Bergeron 1995; 
Saffer, Livio & Yungelson 1998), concluding that about 10% 
of presently known WDs are He WDs. 

We have used our atmosphere code to study the changes 
in colours and absolute visual magnitude of He WDs dur- 
ing the cooling. The calculations were carried out by us- 
ing the cooling sequences of Benvenuto & Althaus (1998) 
for helium-core composition and hydrogen envelopes. The 
chosen evolutionary models cover stellar masses M — 0.15, 
0.30 and 0.45 Mq with Mh/M = 4 x 10"^ 2 x 10"^ and 
4 X 10"*, respectively, where Mh is the mass of the H enve- 
lope. We wish to point out that recent studies (Driebe et al. 
1998; Sarna, Antipova & Muslimov 1998; Sarna, Antipova 
& Ergma 1999) predict the formation of He WDs with mas- 
sive hydrogen envelopes. Since the H-burning increases with 
Mh, thus implying longer evolutionary times, the ages esti- 
mated in this paper (using models with rather thin hydrogen 
envelopes) should be considered as lower limit values. 

We summarized here the most important results. Inter- 
esting observations can be inferred from the analysis of Fig. 



12, which shows evolutionary tracks along with an observa- 
tional sample of cool stars in {J~H, V~K) and {J—K, R—I) 
two-colour diagrams. The cooling tracks in both diagrams 
show an almost linear behaviour following the observed pho- 
tometry up to a predicted effective temperature of approx- 
imately 4500 K. Although there is an important scatter of 
observed colours, it is clear that J — H and J — K show a 
shifting to the blue in the cool extreme of the sample, for 
DA and non-DA stars. This effect is very well represented 
by the theoretical cooling sequences, which gets dramati- 
cally bluer in these colour indices due to the onset of the 
H2 collision-induced absorption for T^n ^ 4500 K. Note 
how the theoretical sequences follow the red edge in V — K 
and R — I oi the observational sample, with the less massive 
models on the right. 

Fig. 12 seems to provide observational evidence of the 
presence of a strong pressure-induced opacity in most of the 
cool atmospheres. It also suggests that very cool WDs iden- 
tified as non-DA, could have an important abundance of 
hydrogen in their atmospheres (perhaps beyond the spec- 
troscopic limit of detection due to the low temperatures) . A 
good example is LHS 1126 (see Fig. 12), the peculiar emer- 
gent flux has been interpreted in terms of H2-He CIA in 
a mixed H/He atmosphere with log 7V(He)/ log 7V(H)= 0.8 
(Bergeron et al. 1994). 

Another feature exhibited in this figure is the predic- 
tion of a second turnover in the J — H colour (around 
V — K — 0.3) for an advanced stage of the cooling. In this 
part of the {J — H,V ~ K) two-colour diagram, the tracks 
of stars with different mass coincide, but the cooling times 
differ significantly. For example, the lowest value of J — // is 
reached by a 0.15 M© He WD at 7 Gyr, while a 0.45 Mq star 
needs 16 Gyr more, due to its lower radius and, therefore, 
to a lower rate of cooling. 

In the {V — I, Mv) colour-absolute magnitude diagram 
of Fig. 13, the tracks corresponding to masses 0.15, 0.30 
and 0.45 Mq are shown from right to left. The diagram 
shows that when a 0.15 Mq He WD is older than 5 Gyr, 
it gets bluer in the V — I colour. There, the absolute mag- 
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nitude changes slowly and the star remains brighter than 
Mv = 16.5 up to 10 Gyr. Similar behaviour is found for 
more massive models, but a longer time (around 20 Gyr at 
M = 0.45 Mq) is spent by these stars before the turnover 
in V — I takes place, as a consequence of an abundant H2 
formation in their atmospheres. 

A number of observed WDs appear redward in V — / 
of the O.45-M0 sequence. Since the stellar evolution the- 
ory predicts that single WDs with very low mass cannot 
be formed within the lifetime of the Galaxy, those stars in 
the diagram correspond to unresolved (or barely resolved) 
double degenerate systems which appear overluminous, or 
single low mass WDs forming a binary system with a much 
less luminous companion. 

We have used our code to derive atmospheric param- 
eters, effective temperature and surface gravity, for five of 
these stars using photometric data of BRL. Combining these 
with evolutionary results we have also derived masses and 
ages for the stars. The effective temperature can be con- 
strained from atmosphere fits to the observed colour indices. 
However, the surface gravity is poorly determined by this 
procedure for hot stars (Te// ^ 4500 K), because the colours 
are rather insensitive to changes of g in such models. Nev- 
ertheless, it is possible to derive the gravity and the mass 
of a star when its absolute magnitude is known. Therefore, 
we determine the effective temperature, the surface gravity 
and the mass of each star, fitting simultaneously its absolute 
visual magnitude and colours with a least-squares method. 
The B — V colour was not included in this procedure because 
there is an opacity source affecting the B filter, which is not 
considered in the model (see Section 4.3). The mass-radius 



relations employed in the fitting and the evolution ages were 
taken from models kindly provided by Benvenuto & Althaus 
(2000, private communication). All interior models used here 
have hydrogen envelopes. 

Treating the observed objects as single WDs, we obtain 
the results summarized in Table ^ Successive columns give 
the WD identifier, name, effective temperature, surface grav- 
ity, estimated mass, mean deviation of the magnitude and 
colour fitting, calculated and observed absolute visual mag- 
nitude, and age. The quoted errors have been determined as 
the changes necessary to produce a mean incertainty of 2a 
provided by the fitting procedure. 

The effective temperatures determined here are on av- 
erage 65 K hotter than those previously estimated by BRL. 
These differences are not significant against to the mean er- 
rors of both calibrations (120 K here and 110 K in the BRL 
work) . As discussed above, the determination of the Tef / is 
mainly provided by the broad-band fit, whereas the temper- 
ature uncertainties can be significantly constrained by the 
absolute magnitude fit. 

Our log p values are 0.1 larger than the results of BRL, 
and we also obtain systematically larger mass for these stars 
(by « 0.06 Mq), in part due to the use of different mass- 
radius relations. Furthermore, using the helium core mod- 
els of Benvenuto & Althaus, we estimated ages about twice 
than the BRL determinations, who used carbon core models 
of Wood^. These differences have an origin mainly in that 



3 Evolutionary models of helium core WDs are available only 
recently. 



the heat content stored per gram in a helium core is larger 
than that of carbon (or carbon-oxygen) core for fixed tem- 
perature. Consequently, He WDs remain bright for compar- 
atively longer times, while the C and C-O WDs cool more 
rapidly. 

Among the five stars analysed here, only L870— 2 is con- 
firmed as binary. This star is a double-lined spectroscopic 
binary composed of a detached pair of DA WDs (Saffer, 
Liebert & Olszewski 1988). LHS 239 is a component of a re- 
solved binary (e.g. Poveda, Herrera, Allen, Cordero & Laval- 
ley 1994). Its companion is LHS 240, a WD classified DC9 
in the McCook & Sion catalog (1987). For LHS 239, BRL 
have found a discrepancy between the effective temperature 
derived from its energy distribution (colour-fit) and that ob- 
tained from spectroscopy (fit to Ha line profile). They con- 
sider that it could be explained if this object is itself an 
unresolved double system. The remain stars in the sample 
(G141 - 2, G187 - 8 and G128 - 7) are suspected double 
degenerate binary (BRL) because their derived masses fall 
below 0.4 Mq and the current theories attribute to such 
WDs a binary origin. 

Table ^ shows the results of interpreting this group of 
stars as unresolved binaries containing two identical WDs. 
The combined mass estimated for these systems is below the 
Chandrasekhar limit, although the individual stars are too 
massive to have helium cores, except perhaps in the case 
of L870 — 2. It is worth noting that combinations of two 
different atmospheres (not analysed here) could allow mixed 
systems of He WD -I- CO WD. The effective temperatures 
in Table ^ are equivalent to those deduced from a single WD 
fitting due to the lack of sensitivity to gravity variations of 
colours, except for LHS 239. 

From a detailed spectroscopic analysis Bergeron, Wese- 
mael, Liebert & Fontaine (1989) have deduced that L870 — 2 
is composed of a pair of DA WDs with comparable effective 
temperatures (Ti = 7470 K and T2 — 6920 K) and masses 
{Ml = 0.47 Mq and M2 = 0.52 Mq). For these determina- 
tions, they adopted the zero-temperature mass-radius rela- 
tion of Hamada-Salpeter (1961). Assuming identical compo- 
nents for the system and using both helium core and carbon- 
oxygen core models recently available, we confirm that their 
individual masses are within the core-helium ignition limit 
(f» 0.45 - 0.50 Mq, see Table |). Our resuh for the effective 
temperature (T^f / = 7200 K) is intermediate between those 
deduced by Bergeron et al.. 

We note that LHS 239 cannot be modelled satisfacto- 
rily by a binary with identical components (the mean de- 
viation of the fit, a = 0.13, is much larger than the un- 
certainty of the observed magnitude, aobs = 0.03), while a 
good fit is obtained considering this star as a single object 
(a — 0.05). These results suggest that the photometry ob- 
served for LHS 239 corresponds to a single He-core WD. We 
consider that another origin, different to unresolved binary 
possibility, would be need to explain the discrepancy in the 
calibrations of T^ff obtained by BRL (see above). 

The photometric fits for G141 - 2 and G187 - 8 do not 
show any clear preference between the single WD or un- 
resolved binary interpretations. In particular, we have not 
found a very satisfactory hydrogen atmosphere fit to the 
observed colour indices for these stars, perhaps due to a sig- 
nificant presence of helium or metals in their atmospheres. 
For G128 — 7, a good fit is obtained for both single and dou- 
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Table 3. Fundamental parameters of low mass white dwarfs assuming a single helium core star 



WD 




Name 


T.ff (K) 


log 9 


M/Mq 


cr (mag) 


cal 


^'W obs 


Age (Gyr) 


0135 - 


052 


L870 - 2 


7200±80 


7.28±0.08 


0.30±0.02 


0.03 


12.41 


12.40±0.07 


1.8±0.3 


0747 + 


073A 


LHS 239 


4210±50 


7.71±0.08 


0.44±0.04 


0.05 


15.65 


15.65±0.03 


11.8±1.5 


1818 + 


126 


G 141 - 2 


6450±200 


7.40±0.11 


0.33±0.03 


0.10 


13.06 


13.04±0.50 


2.2±0.1 


2048 + 


263 


G 187-8 


5250±150 


7.46±0.15 


0.34±0.06 


0.12 


14.14 


14.12±0.15 


4.0±1.3 


2248 + 


293 


G 128 - 7 


5600±80 


7.55±0.07 


0.37±0.03 


0.05 


13.94 


13.94±0.19 


3.3±0.2 



Table 4. Parameters derived assuming a binary with two identical degenerate stars 



WD 


core 


T.ff (K) 


log 9 


M/Mq 


o-(mag) 


My (binary) 


Age (Gyr) 


0135 - 052 


He 


7200±80 


7.81±0.05 


0.50±0.02 


0.03 


12.41 


2.5±0.2 




CO 


7200±80 


7.78±0.04 


0.47±0.02 


0.03 


12.40 


l.liO.l 


0747 + 073A 


CO 


4180±70 


8.22±0.16 


0.72±0.10 


0.13 


15.65 


8.6±1.0 


1818 + 126 


CO 


6450±200 


7.92±0.10 


0.54±0.06 


0.09 


13.04 


1.7±0.4 


2048 + 263 


CO 


5250±170 


7.98±0.15 


0.57±0.09 


0.12 


14.12 


3.8±1.1 


2248 + 293 


CO 


5600±80 


8.09±0.06 


0.64±0.04 


0.05 


13.94 


3.5±0.2 



ble star possibilities. Unfortunatly, for the deduced effective 
temperature (5600 K), the energy distribution is gravity- 
insensitive and the fits do not highlight the binary nature of 
this star. Additional spectroscopic data and line calculations 
could improve the energy distribution fits and the identi- 
fication of the components of these and other (suspected) 
unresolved binaries. 



5 CONCLUSIONS 

We have computed a grid of hydrogen model atmospheres 
for white dwarfs covering an extensive region in the {Teff,g) 
plane. For this, we have employed a partial method of lin- 
earization of equations for stellar atmospheres in which the 
energy is transported by radiation and convection. In view 
of numerical instabilities usually found to solve the resulting 
set of equations, we use a suitable initial temperature dis- 
tribution which improves the convergence properties taking 
into account the low values of temperature gradient expected 
in deep layers of convective models. 

From a relatively simple thermodynamic model, we 
have derived analytical expressions for the adiabatic tem- 
perature gradient and heat capacities of a mixture of hydro- 
gen species. They result appropriate and easy to apply in 
the calculation of the convective flux. In particular, for the 
density regimes present in most of the atmospheres studied, 
our expression of the adiabatic gradient successfully repro- 
duces the values of the ideal equation of state developed by 
Saumon & Chabrier. 

We have analysed in detail the structure and emergent 
flux of hydrogen white dwarf atmospheres over a wide range 
of effective temperatures and gravities. This analysis com- 
plements the most recent studies, to provide an overall pic- 
ture of these atmospheres and the changes to them expected 
during the cooling. In particular, we find strong modifica- 
tions of the emergent flux for Te// ^ 4500 K produced by 
the pressure-induced absorption of molecular hydrogen. The 
results obtained from our code have been compared to those 
recently published. Our calculations confirm the importance 
of infrared photometry for the study of cool white dwarfs. 



and its use in the determination of ages of the Galaxy and 
globular clusters. 

Spectral evolution of helium core white dwarfs has 
briefly been studied using evolutionary sequences recently 
available for these objects. Colour-colour diagrams show 
that the theoretical cooling sequences follow the low mass 
boundary at the cool extreme of an observational sample 
of white dwarfs. It is very interesting to note that on a 
timescale less than 9 Gyr, very low mass He WDs (0.15 
Mq) show strongly non-blackbody colours with absolute vi- 
sual magnitude that are still below 16.5, while massive He 
WDs (« 0.45 Mq) remains brighter than My = 17 for 25 
Gyr. 

Combining atmosphere calculations with interior model 
results, we determined the possible fundamental parameters 
of some low mass white dwarfs. Because, in general, very 
low mass WDs are suspected to be close binaries, the es- 
timated parameters depend on the possible components of 
each system. With the assumption that both components are 
identical, the analysed stars present preferentially carbon- 
oxygen core compositions but with combined masses below 
the Chandrasekhar limit. 

We will present self-consistent cooling calculations in 
the near future, by incorporating our atmosphere code as 
an appropriate surface boundary condition for evolutionary 
models of white dwarfs. 
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• Figure 1.- Comparisons of values for the adiabatic 
temperature gradient Va = 91nT/91nP|g, obtained from 
SCVH EOS (dotted lines) and our expression (Eqs. (|-|), 
solid lines) calculated for different atmosphere models at log 
g = 8 (with effective temperatures indicated on the plot). 

• Figure 2.- The structure of the convection zones for 
models with log g = 6.5 and 8.0. Curves show the layers 
where convection carries 1, 10 (20) 90 and 95 per cent (from 
bottom to top) of the total flux at log g = 6.5 (solid lines) 
and 8.0 (dashed lines). The dotted curves show from left to 
right 95, 50 and 5% of H2 abundance, and 5, 50 and 95% of 
H"*" abundance for models at log g = 6.5. 

• Figure 3.- Temperature-pressure stratifications for at- 
mosphere models with log g = 6.5 (solid lines) and 8.0 (dot- 
ted lines). Models shown have T^// = 2500 (250) 4000, 4500, 
5000 (1000) 12000, 15000 (5000) 30000 (10000) 60000 K 
(from bottom to top). For models at log g — 6.5, the tri- 
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angles indicate the r = 2/3 level, and the circles show the 
regions where convection carries 1% of the total flux near 
the transition between radiative and convective zones. The 
plasma phase transition predicted by Saumon & Chabrier 
(1992) is represented by the heavy curve labeled "PPT". 
This figure extends the results in Fig. 2 of Bergeron et al. 
(1995a). 

• Figure 4.- Runs of temperature with Rosseland mean 
optical depth for models presented in Fig. 3. Solid circles 
indicate the top of the convective zone of models at log g = 
6.5. 

• Figure 5.- Density proflles for models as those referred 
to in Fig. 4. Visible extremes of each curve are located at 
log r = —6 (botton) and 1.5 (top). 

• Figure 6.- Emergent flux distributions for various ef- 
fective temperatures (indicated on the plot) and log g = 6.5 
(solid lines), 8.0 (dashed lines) and 9.0 (dash-dotted lines). 
Those are compared with a blackbody spectrum at the same 
Teff (dotted lines). 

• Figure 7.- In an obvious notation, individual contri- 
butions to the opacity at a depth r = 1 of the models 
with log g = 8.0 and T^// = 60000, 10000, 6000 and 2500 
K, from (a) to (d). These correspond to (log T, log p) = 
(4.76,-6.53), (4.03,-5.75), (3.80,-4.15), (3.53,-1.13), re- 
spectively. Heavy lines in (d) indicate the most important 
opacity contributions at r = 0.1 of the 2500 K model, where 
(logr,logp) = (3.09, -2.46). 

• Figure 8.- Normalized profiles of emergent fluxes for 
Teff = 2500 (dashed line) and 10000 K (solid line) at 
log g = 8, and the blackbody profile at 2500 K (dash- 
dotted line). Dots indicate the filter transmissions of BVRI 
Johnson- Cousins and JHK Johnson- Glass systems. 

• Figure 9.- The (V ~ R,R — I) two-colour diagram 
shows hydrogen model atmospheres of white dwarfs. Heavy 
solid lines represent our models at log g = 6.5 (0.5) 9.0 
from top to bottom, while thin lines connect those with the 
same Te// (labeled). It is also shown the log g — 8.0 se- 
quences of Saumon & Jacobson (1999; triangles) and Berg- 
eron (2000; diamonds); and the cooling track for a O.6-M0 
carbon-oxygen WD calculated by Richer et al. (1999; dashed 
line). Photometric data of DA (filled circles) and non-DA 
(open circles) stars are taken from BRL, except the colour 
indices of LHS 3250 which are taken from Harris et al. 
(1999). 

• Figure 10.- The {B — V, V — K) colour-colour diagram 
shows the photometric data of BRL (circles), the log g = 8.0 
sequence of Bergeron (diamonds) and our model grid (lines). 
Notation is as Fig. 9. 

• Figure 11.- Same as Fig. 10 for the (V - I,V - K) 
diagram. 

• Figure 12.- Evolutionary tracks of He WDs with 
masses 0.15, 0.30 and 0.45 Mq (solid lines from right to 
left) displayed in {J - H,V - K) and [J - K,R- I) two- 
colour diagrams. A subset of cool WDs observed by BRL 
is represented with filled (DA) and open circles (non-DA 
WDs). Ages in Gyr are labeled for 0.15 and O.45-M0 mod- 
els. Dashed lines indicate the cooling track of a black body 
(independently of its size!). 

• Figure 13.- The colour-absolute magnitude diagram 
shows the cooling tracks of WDs with helium cores, hydro- 
gen atmospheres and masses 0.15, 0.30 and 0.45 Mq (solid 
lines from right to left). Ages are indicated in units of 10^ 



years. Dotted lines show the same cooling tracks assuming 
stellar radiation as black body. The observed photometry 
is taken from BRL (DA WDs: filled circles; non-DA WDs: 
open circles); Monet et al. (1992; diamonds) and Harris et 
al. (1999; LHS 3250). 
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